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Abstract 
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Finally a simple model to take into account finite detector resolution is discussed. 
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Whether a source emits photons or massive particles, if it is to be used in an interfero- 
metric experiment, an essential property is its coherence. The study of coherence in optics 
has shown that more than one kind of coherence can be defined l| . The most familiar type 
of coherence is known as first order coherence and is related to the visibility of interference 
fringes in an interferometer. It is proportional to the value of the correlation function of the 
associated field. Second order coherence is less intuitive and corresponds to the correlation 
function of the intensity or squared modulus of the field. From a particle point of view, 
second order coherence is a way of quantifying density correlations and is related to the 
probability of finding one particle at a certain location given that another particle is present 
at some other location. Particle correlations can arise simply from exchange symmetry ef- 
fects and exist even when there is no interaction between the particles. This fact was clearly 
demonstrated in the celebrated Hanbury Brown Twiss experiment which showed a second 
order correlation for photons coming from widely separated points in a thermal source such 
star Q. 

Analogous correlations in massive particles have also been studied, particularly in the 
field of nuclear phvs.cs Q fl B fl Q. Spatial correlations us,ng low energy electrons have 

also been studied [8|, |9|. The advent of laser and evaporative cooling techniques has also 



made it possible to look for correlations between neutral atoms and recently a wide variety of 

different situations have been studied Q Q, Q, Q, Q, H, Q • 

Correlation phenomena are 



generally richer when using massive particles because they can be either Bosons or Fermions, 
they often have a more complex internal structure and a large range of possible interactions 
with each other. In the field of ultra-cold atoms, the many theoretical papers to date have 
included treatments of bosons in a simple three dimensional harmonic trap 



bosonic cloud in the Thomas Fermi regime and Tonks- Girardeau limit 



19|, 



20 



ISLa ID 
2l|, the 

Mott-insulator or superfiuid phase for atoms trapped in optical lattices [22| and the 2D gas 
23|. 

Almost all these theoretical treatments have dealt with atomic clouds at thermal equi- 
librium. On the other hand, all the experiments so far except Ref. 16| have measured 
correlations in clouds released from a trap which expand under the infiuence of gravity 
and possibly interatomic interactions. It is generally not trivial to know how the correla- 
tion properties evolve during expansion. Moreover, matter waves have different dispersion 
characteristics than light. All this raises interesting questions concerning the value of the 
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correlation lengths during the atomic cloud expansion. In particular we would like to know 
how to use the results of Ref.[l3] to analyze the experimental results of Ref. a con- 
ceptually simple experiment in which second order correlations were measured in a freely 
expanding cloud of metastable helium atoms. The correlation length was defined as the 
characteristic length of the normalized second order correlation function. We will use the 
same definition in this paper (see section [TAl for details). 

To illustrate a more general question that comes up in thinking about the coherence 
of de Broglie waves, consider a beam of particles with mean velocity v hitting a detector. 
Two obvious length scales come immediately to mind, the de Broglie wavelength h/{mAv) 
associated with the velocity spread Av and the length associated with the inverse of the 
energy spread of the source hv/m{AvY. These two scales are obviously very different if v is 
large compared to the velocity spread. In this paper, we will show that in an experiment such 



as 



15l | , the correlation length corresponds to neither of the above length scales, although they 



can be relevant in other situations. We find that the correlation length after an expansion 



time t of a cloud of initia 



size s is ht/ms. This result is the atom optical analog o 



Cittert-Zernike theorem 2J]. It has also been stated in a different form in Ref. 251]. For 



the van 



the special case of an ideal gas in a harmonic trap of oscillation frequency uj, the correlation 
length can be recast as Xut where A is the thermal de Broglie wavelength. Hence the 
correlation length after expansion is simply dilated compared to that at equilibrium with 
the same scaling factor as the spatial extent of the cloud itself. 

We will confine ourselves here to the case of a cloud of non interacting atoms released 
suddenly from a harmonic trap. The paper is organized as follows. We will begin in section 
[I with some simple definitions and general results about the correlation properties of a non- 
interacting cloud both at thermal equilibrium in a trapping potential and after a ballistic 
expansion. Without making any assumptions about the form of the trapping potential, 
we can only find simple analytical results in the limit of a non-degenerate gas. Next we 
will make a more exact and careful treatment by specializing to the very important case 
of a harmonic potential. We introduce the flux operator [26] involved in the experimental 
electronic detection with metastable helium and then calculate the correlation function of 
the flux. We will summarize the results and give a physical interpretation in section lllli 
This interpretation will allow us to comment on the rather different case of a continuous 
beam as in the experiments of Ref. [3, [lo|, 14 1. In section HVl we will use our results to analyse 
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the experimentally important problem of finite detector resolution. Finally, the appendix 
adds some detailed calculations concerning the expressions found in section IIIBI 



I. GENERAL RESULTS ON CORRELATION FUNCTIONS OF NON- 
INTERACTING GASES 



Here we recall some basic results concerning the density and first and second order cor- 
relation functions for a cloud of non-interacting bosons at thermal equilibrium. A more 
detailed analysis can be found in Ref.ld]. Theoretical treatments that take into account 



interatomic interactions can be found in Ref. 



17 



27|. We also give some approximate 



results for a non-interacting gas after it has expanded from a trap. 



A. Definitions 



Consider a cloud of atoms at thermal equilibrium at a temperature T, confined in a 
trapping potential. This potential is characterized by {ej,4'j{r)} the energy and wavefunc- 
tion of level j (here supposed non-degenerate for simplicity). In second quantization, one 
defines the field operators 

^t(r) = 5:^f(r)aj, ^{r) = J2Ur)ay 
j j 

The operator aj creates and dj annihilates one particle in state whereas ^"''(r) creates 
and ^(r) annihilates a particle at position r. 

Correlation functions and the atomic density are statistical averages of such field opera- 
tors. We use the Bose-Einstein distribution, (ajak) = Sjk{e^^''^^^^ — 1)^^ where /3 = l/i^ksT), 
ks is the Boltzmann constant and fi is the chemical potential. The value of /j, ensures the 

normalization J2{0']<^j) = We can then define 
j 

• the first order correlation function G^^\r,r') = (^"''(r)^(r')), 

• the second order correlation function G(^)(r,r') = (^'''(r)^'(r)^'^'(r')^'(r')) 

• and the density Peg(r) = (^l'(r)^(r)) = G(^)(r,r). 
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Several other first and second order correlation functions can be defined (see below) but 
these are the most common ones. The first order correlation function appears in interference 
experiments whereas second order correlation functions are related to intensity interference 
or density fluctuation. First and second-order correlation functions are connected for thermal 
non-interacting atomic clouds. The G^"^^ function contains a statistical average of the type 
^a\^a\an) which can be calculated through the thermal averaging procedure (Wick theorem 



28|). One finds (ajakciian) = (ajaj)(akak)(5ji5kn + ^jn^ki) + (ajaj)5ki5jn, which leads to 

G(2)(r, r') = pe,{v)pe,{v') + |G(i)(r, v')\^ + Pe,(r)5(r - r') 

The last term is the so-called shot-noise term. It will be neglected in the following because 
it is proportional to whereas the others are proportional to A^^. 

It is convenient to define a normalized second order correlation function 

(2), n G(2)(r,r') 



Peg(r)peg(r')' 

If the cloud has a finite correlation length, then for distances larger than this length the first- 
order correlation function vanishes. Then (7*^^^(r, r) = 2 and (^'^^^(r, r') 1 when |r — r'| — »• 
oo. This means that the probability of finding two particles close to each other is enhanced 
by a factor of 2, compared to the situation where they are far apart. This is the famous 



bunching effect first observed by Hanbury Brown and Twiss with light [2|. 

The above expression of the G^^-* function cannot be applied in the vicinity and below 
the Bose-Einstein transition temperature. The calculation of (aja-kQ-ifln) is performed in the 
grand canonical ensemble which assumes the existence of a particle reservoir that does not 



exist for the condensate. It is well known 



29| that this gives unphysically large fluctuations 



of the condensate at low enough temperature. This pathology disappears at the thermody- 



291]. It has also been shown that it cancels 



namic limit if there is an interatomic interaction 
for a finite number of non-interacting particles if one uses the more realistic canonical en- 
semble 30|. One way to keep using the grand canonical ensemble is to add the canonical 



result for the ground-state [171]. This approach is validated by the results in Ref. SOj and 
will be used in the following. The largest deviation is expected to occur near the transition 
temperature [s^. The contribution of the ground state is — (aoao)^5jo(5ko'^io'^no- Then, with 
Po the ground-state density, it follows that, 
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G(2)(r,r') = pe,(r)pe,(r') + |G«(r, r')!^ - Po(r)po(r') (1) 
The normalized second order then becomes 

Peg(r)Peg(r') Peg(r)peg(r') 

Because the ground state density is neghgible for a thermal cloud, the normalized cor- 
relation function g^'^\r,r') still goes from 2 to 1 as the separation of r and r' increases. 
On the other hand, for a BEG at T = 0, only the ground-state is occupied. Then 
\G^^\r,r')\'^ = Peg (r) Peg (r') = po(r)po(r') and g^'^\v,r') = 1. The amount of particle bunch- 
ing present in the second order correlation function can be quantified as g^'^\r,r') — 1 and 
this typically decays exponentially as the modulus squared of the separation between the 
two points increases. We define the correlation length to be the characteristic length over 
which the amount of particle bunching decays, that is the distance over which g^'^^v, r') — 1 
decays to 1/e of its maximum value. The correlation length of a BEG is infinite. Such a 
system is said to exhibit bunching at high temperature over the correlation length and no 
bunching in the condensed phase. 



B. Correlations in an expanding cloud 

In most experiments, particle correlations and other characteristics are not directly mea- 



sured in the atom cloud, (Ref. [l6| is an exception). Rather, the cloud is released from a 
trap and allowed to expand during a "time of fiight" before detection. For a sufficiently long 
time of fiight, and neglecting interactions between the atoms, the positions one measures at 
a detector refiect the initial momenta of the particles. The results of section II Al concerning 
the correlation functions in position space all have analogs in momentum space. In fact the 
correlation functions in the two reciprocal spaces are closely related. At equilibrium, i.e. 
inside the trap, the following relationships can be easily derived: 

I dp G^^\p, p)e-'P"/^ = JdK G^^\K - r/2, R + r/2) 

J dr (r, r)e^'i-"/^ = / '^^^^ " ^/^' ^ + ^/^) 
In other words, the spatial correlation length is related to the width of the momentum 
distribution and the momentum correlation length is related to the width of the spatial 



distribution i.e. the size of the cloud. No equally simple and general relationship holds 
for the second order correlation functions. This is because, close to the BEC transition 
temperature, and at points where the ground state wave function is not negligible, the 
special contribution of the ground state, the last term in Eq. [T] must be included, and this 
contribution depends on the details of the confining potential. On the other hand, for an ideal 
gas far from the transition temperature one can neglect the ground state density, make the 
approximation that the correlation length is very short, neglect commutators such as [f, p], 
and then write the thermal density operator as 6" = e~^2^e~^^^^\ These approximations 
lead to: 

G(^nP,P') =Pe.(p)Pe,(p') + IG^^HP-POP 

and, 

G^^\P - q/2, P + q/2) ~ e~^^ J dr e-^^^^'^e'^ 

One sees that in this limit, the interesting part of G*^^^ in momentum space is proportional 
to the square of the Fourier transform of the density distribution and independent of the 
mean momentum P. This result is the analog of the van Cittert-Zernike theorem {2^. For a 
trapped cloud of size Sa in the a direction, one has a momentum correlation "length" given 
by: 

= -• (2) 

If atoms are suddenly released from a trap and allowed to freely evolve for a sufficiently long 
time t, the positions of the particles reflect their initial momenta and the spatial correlation 
length at a detector is given by 

l^^ = = (3) 

m msa 

The normalized second order correlation function is then a Gaussian of rms width l^'^^ / \pl. 
This result was experimentally confirmed in Ref . [l^ . One wonders however, to what extent 
the approximations we have made are valid. The clouds used in Ref. ISj were in fact very 
close to the transition temperature so that effects due to the Bose nature of the density 
matrix may be important. Although the time of flight was very long, it is useful to quantify 
the extent to which identifying the momentum correlation length in the trap with the spatial 
correlation length at the detector is accurate. Finally, the effect of gravity on the falling 
atoms never appears in the above approximate treatment, and we would like to clarify the 
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role it plays. In order to answer these questions we undertake a more careful calculation. We 
will confine ourselves to atoms initially confined in a harmonic trap, a good approximation 
to the potential used in most experiments, and happily, one for which the eigenstates and 
energies are known exactly. 



II. DENSITY AND CORRELATION FUNCTIONS FOR A HARMONIC TRAP 



A. At equilibrium in the trap 

The eigenf unctions for a 3-dimensional harmonic potential of oscillation frequency uja in 
the a direction, are given by: 



a=x,y,z 



Here (Tq = y;^^ is the harmonic oscillator ground-state size, Hj^ is the Hermite poly- 
nomial of order ja and Aj^ = ( x/TrcTrv 2-^" ( ?n- ) ~ • The eigenenergies are given by 



J2 h^aUa + 1/2)- Then [l7|, l29|], with = I3huja and /i = /x — TiY.^al'^, one 

a=x,y,z 



finds: 
and 

G'W(r,r') 



TT 



-| oo 



■tanh(Ifi) 



3/2 



TT 



1 oo 



/=1 



1 



3/2 



exp 



The above expressions can be transformed into more familiar forms in limiting cases: 



For high temperature, /x — — oo and one recovers the Maxwell-Boltzmann distribu- 

— 1^ Q I 

tion. The density is Peqi^^) = ^YlTa^ ^ with A = J^^^^j. the thermal de Broglie 
wavelength. The size of the cloud is Sa = c^al \ff^ - 

The first order correlation function is 



kBT_ 



N 



(4) 



Using our definition, the correlation length is Z*^*^ = 



TT. 
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For a temperature close to but above the Bose- Einstein transition temperature, one has 

— ^ 

to keep the summation over the index /. The density is Peg(r) = p'fi'3/2[e^'^ 11 e ^ "^"j 

a 



oo 



where ga{x) = 2; //'^ is a Bose function. The first order correlation function is 
1=1 

/\ 1 ^e^'^f'f-f _IaL(^Si±lk)2 _^IL(IR^)2 

As the temperature decreases, the number of values of / that contribute significantly 
to the sum increases. It is then clear from the above expression for G^^'^ that the 
correlation length near the center of the trap will increase and that the normalized 
correlation function is no longer Gaussian. Far from the center, only the / = 1 term is 
important and the correlation function remains Gaussian. Thus close to degeneracy 
the correlation length is position-dependent (for an explicit example see Sec JII B 5|) . 



Near and below the transition temperature, the second order correlation function 
is given by Eq. ([1]) with po(r) = jzr^ U f^^-^ ■ As the temperature decreases, 
the correlation at zero distance, g^'^\0,0) decreases from 2 to 1 and the correlation 
length increases. Around the transition temperature, g^'^\0,0) is already significantly 
different from 2 since the condensate peak density is already very large for a non- 



interacting harmonically trapped cloud 
infinite and g^'^\r,r') = 1. 



31| . At T = 0, the correlation length is 



B. Correlations in a harmonically trapped cloud after expansion 

Here we consider the cloud after expansion. First we discuss two classes of detection 
methods which must be distinguished before calculating correlation functions. 



1. Detection 

We assume that the trapping potential is switched off instantaneously at t = 0. The 
cloud expands and falls due to gravity. Two types of detection can be performed: 

• Snap shot. An image is taken of the entire cloud at t = to- We have then access to 
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The usual imaging technique is absorption, and so one has access to the above cor- 
relation functions integrated along the imaging beam axis. This was used for the 
experiments of Ref. 



12 



13|. 



Flux measurement. The atoms are detected when they cross a given plane. We will 
only consider the situation in which this plane is horizontal at z = H. One has access 
to 

Gfl (r = {x, y, z = H}, t; r' = {x', y', z' = H}, f) = (/(r, t)/(r', t')) 

where / is the flux operator defined below. The detectio n sy stems required for such 



experiments correspond most closely to those of Refs. 



10 



15| . in which a micro- 



channel plate, situated below the trapped cloud, recorded the arrival times and in 
one case the positions of the atoms. It also corresponds closely to imaging a cloud 
that crosses a thin sheet of light 32 1, or to the experiment of Ref.|l^, in which the 



transmission of a high finesse optical cavity records atoms as they cross the beam. 

These two correlation functions are different, but if the detection is performed after a long 
time of flight, they are in fact nearly equivalent. This equivalence will be discussed in the 
following. 

The flux operator is defined quantum-mechanically by 



J(r,t) = —Im\¥{Y,t)d^i>{r,t) 
m L 



2mi 



¥(v,t)d,i>{v,t)-d,¥{v,t)i>{v,t) 



The flux has thus the dimensions of a density times a velocity. We will give the explicit 
expression of this velocity in the section fll B 4[ Here, the atomic field operators ^(r, t) 
depend on space coordinates as well as on time. They represent the time evolution of the 
atomic field during the flight of the atoms, falling from the trap. The field operators for the 
falling cloud can be easily derived if we assume that there are no interactions between the 
atoms and that the occupation number in each mode is constant (as in free expansion). In 
this case, these operators can be defined as 

3 j 

where the spatiotemporal dependence is carried by the wave function and the statistical 
occupation by the creation and annihilation operators. 
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2. Ballistic expansion of a harmonic oscillator stationary state 



After switching off the trap, the harmonic oscillator wave-functions noted ip^ are no 
longer stationary states. There are two ways to calculate the correlation after expansion: 
propagation of wavefunctions or propagation of the density matrix (the Schrodinger or the 
Heisenberg picture). In the following we will use the first approach which is physically more 



transparent (see 



33j for the Heisenberg picture) 



The ballistic expansion of a cloud is easy to calculate with the appropriate Green function. 
The Green function K is defined as 



/oo 
dro K{r,t;ro,to) V'j'(ro,to)- 
-oo 



As the V'? functions are stationary states for t < 0, we can take to = in the following. The 



Green function for particles in an arbitrarily time- varying quadratic potential is known 34 1. 
After expansion, the potential is only due to gravity and the Green function is then 



( ^ \ Ja[T-ro)'^ ih{z-^Zi)) -ic 



witha = i^,6=^andc 



24h ■ 



One can then derive an analytical expression of ^j(r, t) 35|, |36|: 



V'j(r,t) 



,i0(r,t) 



n 



„ija{Sa+3n/2) 

.0 



^/uJat 



where 5a = tan ^[^1, 
and, with f = {x, y, z}, 



[xu;.jf + {yuoytf + {iuo.tf + 2ge{z - -gf 



X 



X 



y 



H - \ge 



377 
T 



(5) 



(6) 



(7) 



'I + ujlt^ + uolt^ ^1 + uolt^ 

In the case of flux measurement, the position of the detector is fixed aX z = H . The 
phase (f){x,y,t) is global as it does not depend on the index j; it will cancel in second 
order correlation measurements. This is in contrast to interferometric measurements where 
it is this phase that gives rise to fringes. The above results show that after release, the 
wavefunction is identical to that in the trap except for a phase factor and a scaling factor 
in the positions [3^. This scaling is obviously a property of a harmonic potential, and it 
considerably simplifies the expression of the correlation functions as we will see below. 
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3. Flux operator 



Using dzHn{z) = 2nHn^i{z), the spatial derivative of the wavefunction can be written: 



where the velocities Vi V2 and are time dependent and are given by 



fl(t) = CU; 



^ 1 + uoH^ 



H + -of 



H - \gf 



1 + C^2t2 



V2( 



(8) 
(9) 
(10) 



'1 + uolt^ 

The velocity f 2 is usually much larger than the other two and will give the dominant con- 
tribution for the mean flux and the second order correlation function. An atom with zero 
initial velocity will acquire after a time t a velocity gt which is close to V2{t). The flux 
operator is, 



/(r,t) = E 



(11) 



where j — Iz is the vector {jx,jy,jz — 1) and where we write = ip{r, t). 



4- Mean density and mean flux 



We will first calculate the mean density p{r,t) = (^"''(r, t)^(r, t)). Using Eq.dS]), one 



finds easily that p(r, t) 



Peqi^)- This means that the density has the same form 



during expansion up to an anisotropic scale factor given by Eq.d?]) |37|, l38|. The statistical 
average of Eq. (ITT]) leads to 



(/(r,t))=E 



V2\lpi 



Because v^ip* ipj 



_M ,1,0 {z\j,0 



Then, without any approximation, 



'\lr,J^z)'>\); _^{z) = —vl'ipj^i']^-!^ the second term cancels out. 



(/(r,t)) 



V2{t) 



Peq{r) = V2{t)pir,t) 
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The flux is proportional to the density of a cloud at thermal equilibrium with rescaled 
coordinates. This means that the mean flux of an expanding non-interacting cloud is pro- 
portional to the atomic density without any approximation. This results holds with and 
without gravity taken into account. 

5. Second order correlation 

Here we calculate the correlation functions. A discussion is given in the next section. 
The snap-shot correlation function is 

j,k,l,n 

Using Eq.([n]), one flnds, without any approximation (except the neglect of the shot-noise 
term) : 

G'S(r,t;r',t) = ^ ^^^^^ (pe.(f)Pe,(fO + ,?')!'- Po(f)po(f)) • 

As in the case of the mean density, the snap-shot correlation function has the same form as 
in the trap except for an anisotropic scale factor. 
The calculation of G'^JI is similar: 

(/(r,t)/(r',t')) = - E [^K^.^k) - X W:{d.^'^) - (S.V^DV^;] (ajakata„) 

Two major differences appear compared to the mean flux calculation: the terms in fs and 
the phase factor 5q, + 37r/2 in Eq.([5]) do not cancel. This makes the exact calculation very 
tedious. It is postponed to the appendix. 

Experiments are usually performed in situations satisfying two conditions: (1) the width 
of the cloud after expansion is much larger than that of the trapped cloud, and (2) the mean 
velocity acquired during free fall is much larger than the velocity spread of the trapped 
cloud. The flrst condition means that oJoX ^ 1 and the second one that gt ^ y/ k^T jm. 
The latter condition also means that the mean arrival time, to = \J'^H/ g, is much larger 
than the time width ^Jk^rjrmp' of the expanding cloud. With these approximations the 
scale factors become quite simple, x ~ , "V ~ and z ~ ^ ~ 9{to-t) ^ particular, 
the coordinate z is proportional to the arrival time t. This means that in experiments that 
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measure arrival times, the results have the same form when expressed as a function of vertical 
position. 



In the correlation function of the flux, the above approximations also lead to f 2 ~ V^^H 



and Ia/tJ f 3/^2 1 ~ \J^J~'^2H ~ '^tn "^^ere Sz is the width of the cloud inside the trap and 
where the typical value of the occupied trap level, jz, is ~ The term containing ^3 is 
then very small compared to the one proportional V2. In Ref. jisl for instance the above 
ratio is ~ 10~^. We will neglect terms containing ^3 in the following. The phase factors 5a 
in Eq.Q are also very small since cj^t ^ 1 and can be neglected (see appendix IVI D 20 . 



Under all these approximations, one finds 

G^?(r,t;r',0 = , ""f'l^ ^ (Pe,(f)pe.(f') + |G«(?,f')P -Po(f)Po(f' 



We again find the same correlation function as in the trap, rescaled by a slightly different 
factor compared to G\„^ . This factor simply reflects the expansion of the cloud between the 
times t and t'. 

The scaling laws for the harmonic potential result in a very simple expression for the 
correlation lengths at the detector: 



4'^) = X v'l + M'- (12) 

Where l'^^ is the correlation length along the a direction at the detector and l^^ is the 
correlation length in the trap. If the gas is far from degeneracy Z^*-* = and we recover 
the result of Eq. [3l Close to degeneracy the correlation length is position dependent. In the 



case of a pulse of atoms as in Ref. 15|], this formula applies along all three space axes. In 
addition, when making a flux measurement, one often expresses the longitudinal correlation 
length as a correlation time. For a pulse of atoms from a harmonic trap, with a mean 
velocity v at the detector, the correlation time is: 

^icoh) ^ ^ lit) ^ /^3) 
V g 

It is independent of the propagation time as long as Uzt 3> 1. 

These calculations are illustrated in the following figures. For simplicity we have used an 
isotropic trapping potential. As pointed out above, the normalized second-order correlation 
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FIG. 1: (Color online) Two-body normalized correlation function at the trap center, g^'^\f,0) 
for 10^ atoms as function of the position r = r/ujt for various temperatures around transition 
temperature. The horizontal axis is labelled in units of the size of the harmonic oscillator wave 
function a. The thick dashed line corresponds to the transition temperature T* defined in Ref. 31 1 



and is 93.37 huj/kB for 10^ atoms. The temperature step is 0.4 Uuj/kB- The thermal de Broglie 
wavelength is ~ 0.26 a. The effect of the ground state population is clearly visible in the reduction 
of (^^^^(0,0), and in the rapid flattening out of the correlation function slightly below T* . 



functions g\J^ and g^l are virtually identical with typical parameters (see appendix IVI Cp 



and we will use the shorter notation g'^'^\ In Fig. [T] we show the normalized correlation 
function 5f*^^)(f, 0) as a function of f ~ r/cut for various temperatures in the vicinity the 
Bose-Einstein phase transition T*. We use the saturation of the excited state population 
to define T* [31]. This is the correlation function at the center of the cloud. One sees that 
at T = T* (the thick dashed line in the figure), the correlation function at zero distance is 
already significantly diminished compared to higher temperatures. The correlation length, 
on the other hand, is larger than Xujt /\/2ti. Also, one sees that the correlation function is 
almost flat for temperatures a few percent below T* . 

In many experiments of course, one does not measure the local correlation function, 
but the correlation function averaged over all points in the sample |l5^. The effect of this 
averaging is shown in Fig. d We plot g^^\f) = j Jw^ff^^^-^cw (r^r) where the vector 
e is a unit vector in some direction. One sees that the amplitude of the correlation function 



15 




FIG. 2: (Color online) Two-body normalized correlation function Qm {f) for 10^ atoms as a function 
of r. This function is an average of the two-body correlation function over the cloud. The conditions 
are the same conditions as for Fig(TJ Unlike FiglH the shape is always almost Gaussian and 
converges more slowly to a flat correlation for low temperatures. This is because only a small 
region around f = is fully sensitive to the quantum atomic distribution. 

decreases more slowly, and that after averaging, the correlation length hardly varies as one 
passes T*. 

To illustrate how local the effects which distinguish Figs. [1] and [2] are, we also plot in 
Fig. [3] the value of g^'^\f,f), the zero distance correlation function as a function of f in 
the vicinity of the cloud center. One sees that even below T*, the correlator is close to 2 
at a rescaled distance of a few times the harmonic oscillator length scale. We can simply 
interpret this effect by observing that at f the effective chemical potential is fi — V{r). Away 
from the center, the effective chemical potential is small and this part of the cloud can be 
described as a Boltzmann cloud. 

Before interpreting these results further, we recall some of our assumptions and their 
possible violation. First, we obtain Eg. (1121) if we make a semi classical approximation 
assuming that ksT greatly exceeds the energy spacing in the trap in each dimension of 
space. In an anisotropic trap, this condition can be violated in one or two dimensions and 
then correlation length along these directions will be larger and can become infinite for a 
small enough temperature. Second, we have assumed a non-interacting gas throughout. 
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FIG. 3: (Color online) Two-body normalized correlation function g^'^\f, f) for 10^ atoms as function 
of f. The conditions are the same as for Fig(TJ Even for T < T* the correlation goes to 2 far from 
the center. This is due to the finite spatial extent of the condensate. It can also be understood 
in terms of the chemical potential /x(f) which, in a local density approximation, decreases as r 
increases and thus the correlation is equivalent to that of a hotter cloud. 

Repulsive interactions inflate the trapped cloud, and thus reduce the length l^'^^ at the 
detector. We expect this to be the main effect for atomic clouds above the Bose-Einstein 
transition threshold, where the effects of atomic interactions are typically small. The reduc- 
tion is typically a few percent. Even slightly below T*, the condensate density is quite high, 
expelling the thermal atoms from the center of the trap. The effects of interactions inside 
the trap and during the cloud's expansion cannot be neglected. Taking them into account 
is then complex and beyond the scope of this paper. 

III. PHYSICAL INTERPRETATIONS 

The main result of this paper is that in an experiment which averages over a detector 
in the sense of Fig. [2], even at T = T*, the correlation lengths at the detector are well 
approximated by: 

= X 
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The correlation length increases linearly with the time of flight. A simple way to understand 
this result is to consider the analogy with optical speckle. Increasing the time of flight 
corresponds to increasing the propagation distance to the observation plane in the optical 
analog. The speckle size, i.e. the correlation length, obviously increases linearly with the 
propagation distance. Another way to understand the time dependence is to remark that 
after release, the atomic cloud is free and the phase space density should be constant. Since 
the density decreases with time as WioJat) and the spread of the velocity distribution is 
constant, the correlation volume must increase by the same factor 25|. 

Yet another way to look at the correlation length is to observe that, far from degeneracy, 
the correlation length inside the trap is the thermal de Broglie wavelength, that is, -4= 



27r 

Tij iS.'p where Ap = mAf is the momentum width of the cloud. By analogy, after expansion, 
the correlation length is h/{Ap)ioc.i where (Ap)ioc is the "local" width of the momentum 
distribution. As the pulse of atoms propagates, fast and slow atoms separate, so that at a 
given point in space the width in momentum is reduced by a factor 

For a continuous beam, the formula (|T2ll only applies in the transverse directions. In 
the longitudinal direction, an argument in terms of a local thermal de Broglie wavelength 
can be used to flnd the coherence length or time. If the atoms travel at velocity v without 
acceleration, the momentum spread and correlation length remain constant. Deflning the 
energy width of the beam as Ai? = mvAv, one flnds a correlation time X/v = h/AE [7|. In 
the presence of an acceleration such as gravity, the momentum spread of the beam decreases 
(the energy spread at any point AE is constant), which increases the correlation length. 
The correlation time, however, remains h/AE 10 1. 

The result that the coherence length of a cloud of atoms can vary with the distance of 



propagation, is in apparent contradiction with the results of Refs. [39|, |40|] . Those papers give 
convincing reasons, both experimental and theoretical, for why the dispersion associated with 
the propagation of massive particles should not result in an increase of the coherence length. 
The contradiction is resolved by noting that the Mach-Zender interferometer considered in 
that work is sensitive to the function /(r, t) = J dH G^^^H, t; R + r, t). If the Hamiltonian 
commutes with the momentum operator, i.e. if plane waves are stationary states, one can 
easily demonstrate that the function / and hence its width are independent of the time t. 
The experiments we analyze are sensitive to the modulus of G^^^ whose width will always 
increase with time. Thus the coherence length can depend on the interferometer as well as 
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the source. 

The role of the acceleration of gravity in these experiments is minor. It governs the prop- 
agation time and the speed of the particles when they reach the detector. In a pulsed beam, 
gravity has no effect on the correlation length, although it does affect the correlation time. 
It also renders the rescaling of the z coordinate linear for large times so that the correlation 
function in position z and time have the same form. Without gravity (cancellation with a 
magnetic field gradient for example), a pulse of atoms would take longer to reach the detec- 
tor, thereby giving the correlation length more time to dilate, and in addition they would 
hit the detector at a lower velocity. The correlation time would then increase with time and 
its order of magnitude would be = -^^to where vt = is the thermal velocity and 

to = vt/H is the time of flight to the detector. 



IV. EFFECT OF FINITE DETECTOR RESOLUTION 



In the preceding sections, the detector was considered ideal, i.e. with arbitrarily good 
spatial and temporal resolution. Here we will consider a model of a more realistic detector, 
in which we suppose that the spatial resolution in the x — y plane is Gaussian. This is often 



the case due to smearing in pixels [13|, ll6|] and is also approximately true in Ref . [151] . To 
simplify the discussion we will restrict our analysis to the case T ^ T* and use a Maxwell- 
Boltzmann distribution rather than Bose-Einstein distribution. In this case, each direction 
of space is independent and we will only consider one direction at a time in the following. 
There are three different scales in the problem: the size of the cloud at the detector 



s{t) ^ V^^^' correlation length at the detector Z^'^) and the r.m.s. width of the 
detector resolution function d. The definition of the resolution function is that for a density 



p(x) = Ae 2l(tj2 ^ observed density is given by a convolution: 

pohs{x) = / dxopixo) 7=— = , =e ^I^w^. 



(x + x')^ (x-x'f 

Similarly if G^^\x,x') = Ae^'^e ^(2sP~e acIMya" jg ^^le first order correlation function and 
the observed one, we have 



lf X-XQ .2 _!(■" 



= dxodx',\G('\xo,x',)\' '' (14) 
J y/zTcd v27ra 
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/,2 



^(i + dVs2(t))(i + 4rfV(/w)2; 

Consequently, with a = x,y and z: 



g 4(s(t)2+d2]g (i(d))2+4d2 (^;]_5^ 



The amphtude of the normahzed correlation function becomes 



^ob.(0, 0) - ( j - 1 + n Y i+4rf2/(4d))2 • 



The observed widths of the cloud are Sait) ^ s'^{t) + da- 

The observed correlation lengths are l^^^ ^/{¥yTi2daf. The factor 2 can be 
understood as x where the first term comes from the fact that da is defined 
for one particle and not for a pair of particles and the second one comes from the fact 
that the correlation length is not defined as an r.m.s. width. 



In the experiment of Ref. 



15| the trapped cloud had a cigar shape. At the detector the 



cloud was spherical but the correlation volume was anisotropic with l^^ <^ d ^ ^'^V^- In 



the third (vertical) direction, the resolution width was much smaller than any other length 

r(d) 

'_X 

2d • 



scale. The observed contrast of the correlation function was therefore approximately, -fj 



V. CONCLUSION 



The most important conclusion of this paper is that the expansion of a non-interacting 
cloud from a harmonic trap in thermal equilibrium, admits a rather simple, analytical treat- 
ment of the time variation of the density and the correlation functions. In such a pulse 
of atoms, correlation lengths scale in the same way as the size of the density profile. The 
agreement with experiment indicates that the neglect of interactions is a good approxima- 
tion above the BEC transition temperature. An important next step however, is to examine 
interaction effects so that the next generation of experiments, which will be more precise 
and better resolved, can be fully interpreted. 
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VI. APPENDIX 



A. Explicit expression of the flux correlation function 



We found in section IIIBt the following expression for the flux operator: 



j,k 

where j — Iz is the vector {jx,jy,jz — 1) and where we write = 



The second order correlation function for the flux is then, 
(/(r,t)/V,t')) = 



j,k,l,n 



X 



Neglecting the shot-noise and ground-state contributions, this leads to 

(/(r,t)/(r',f)) = (/(r,t))(/(r',t')) + Re(A) 
with 



We write A = J^Ti where the Tj terms can be recast, using tan6a = 1/uJat, tan 5^ = l/uat', 
A, = 5; - S^EjaiS'^ - 5a) = j.A, = ^O(f) and = ^P^{^'). 

a 



Ti = V2V2 E V'j i^l^j) i^l^i) 



a 

'"2V^ 



T2 = l^^st-sE v^7XV^jV^j_i^^i-i>'i(«j«j)(«iai) 
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^3 = Ivsv'l J2lz^ifj^i-u^'*i-u{alaj){alai) 



a 

j.i 

a 

• T5 = -v'2V3 E V^V'j'^jVl-lzV^'l(«loj)(«lOl) 

= -' iw..4^U^i.-, (?>*>'J^-"(aJaj)) (EvC<-,#fe— (a!a,> 

a 

The term Ti is a real number which is not the case for T3, T4 and T5. 



B. Calculation for harmonic oscillator stationary states 

All the above terms can be calculated analytically. All the series are identical in the 

direction x and y. We are then left with the calculation of three series in only one direction: 

00 

n=0 

00 

n=0 

00 

n=0 



The function gu{z, z') = E 4'ni^)'^ni^')^ is known [l7|,l29|] and its expression is gu{z, z') 

n=0 

^exp[-tanh(f) (^)' - coth(f) 



(T^7r(l-e-2") 

Using ~zi)l{z) = -^{z\a + at|^0) = ^[^^o_^(5) + + one finds 

It follows easily that 



00 



n=Q 

00 



n=0 

00 

Moreover, E n^^_^{z)tP^_^{z')e-''^ = e-^[gu{S,z') - dugu{z,z')]. Then, 



71 = 

71=0 2 2 
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C. Explicit expression of the flux correlation function-Part II 



We define Gy(r, r',u) = J2'^ni^)'4^n{^')c~'^"- This function, the 3D equivalent of 

n 

the function g^, is connected to the one-body correlation function by G'^^-'(r,r') = 
E e^'^CgV, r', It) with r„ = /J/ic^^. 



1=1 



Then, 



Ti 



To 



X 



a 

2^^ <j l_e2('=^2+»A^) '-^B 11-7 I , ft-T -|- 



|d3w' I 



T4 = —2 



V2\v'\ 



X 



a 



Ee^'^^|ES^G«(f,F,/T-zA; 



X 



The dominant term is Ti and is the one used in section III B 51 



D. Contribution of neglected terms in the correlation of the flux 

Here we evaluate the neglected the terms T2 to T5 and the shot-noise contribution. They 
will be evaluated in the case of clouds far above BEC threshold. Under this assumption, all 
the functions are separable in the variables x, y and t and the summation over the index / 
in the previous equations reduces to the single term / = 1. 



23 



1. Shot-noise contribution 

Using the above analysis one can show that the main term is still proportional to V2V2. 
The additional term is then, 



V2V2 



=e^'^Gi^)(r,?',r-iA)Gg)(f,f',zA) 



For t = t', A = and G^''(r, r', 0) = 6(y — r'). The shot-noise term is then 

a 

As expected, this term corresponds also to the one at equilibrium with rescaled coordi- 
nates. 

2. T2 — T5 contribution 

We have G}.?(r, t; r', t') = (J(r, t)J(r', t')) = (/(r, t))(J(r', t')) + Re(A) where A = E 
• Case t = t' . 



A = 0, 

2 



then Ti = ^ , "'"^ 
Re(T4) = Re(T5) = 0. 



, T2 and T3 are real number and 



One finds, to leading orders, 

^(2)(o,0,t;0,0,t) - 2 ^ I (f)^(l -2^)(1 - ^) where is the initial size of 



the cloud in the vertical direction and to = y2H/g. 

The deviation from 2 is extremely small in the experimental conditions of [l^ 
(~ 10^^^) but shows that the bunching is strictly speaking not 2 at the center. 
This behavior is expected for any fiux correlation function of dispersive waves 



4l|. 



— The correlation lengths at the detector are not modified by the additional terms. 
Case t^t'. 
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The correlation function can be written as 




and 



^-I(t)'[l-(^)](1-|) 



3 



We have neglected terms in r^, (t-to)^ (^'-^o)^ {t - to^ {f - h) , {t -to){t' -tof 
and higher orders. 

— The value of e is extremely small (~ 10~^°) using Ref . [isl . The deviation from 



the value of the correlation time in the wings of the time-of- flight by ~ 3 %. The 
effect of the phase factor A is negligible. 
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